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Abstract 

Using molecular dynamics simulations, we investigate the relation between the dynamic tran- 
sitions of biomolecules (lysozyme and DNA) and the dynamic and thermodynamic properties of 
hydration water. We find that the dynamic transition of the macromolecules, sometimes called 
a "protein glass transition", occurs at the temperature of dynamic crossover in the diffusivity of 
hydration water, and also coincides with the maxima of the isobaric specific heat Cp and the tem- 
perature derivative of the orientational order parameter. We relate these findings to the hypothesis 
of a liquid-liquid critical point in water. Our simulations are consistent with the possibility that 
the protein glass transition results from crossing the Widom line, which is defined as the locus of 
correlation length maxima emanating from the hypothesized second critical point of water. 

PACS numbers: 61.20.Ja, 61.20.Gy 
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Both experiments and computer simulation studies have shown that hydrated proteins 



undergo a "glass-hke" transition near 200 K 



IBS 



above which proteins exhibit dif- 



fusive motion, and below which the proteins are trapped in harmonic modes. An irnportant 



mum. 



issue is to determine the effects of hydration water on this dynamical transition 
Experiments and computer simulations suggested that when a protein is solvated, the pro- 
tein glass transition is strongly coupled to the solvent, leading to the question of whether 
the protein glass transition is directly related to a dynamic transition in the surrounding 
solvent 111 . 

It has been hypothesized that the deeply supercooled region of the phase diagram of 
liquid water may contain a first order phase transition between a high density liquid (HDL) 
and a low density liquid (LDL) This line of phase transition has a neg ative slope in 



the P-T phase diagram and terminates with a critical point, jl^ 13 . Q, Q| 



•rf 



which is 



located at Tc2 = 200 ± 20 K and Pc2 = 250 ± 100 MPa fig. Upon crossing the first 
order phase transition line above the critical pressure, the thermodynamic state functions 
change discontinuously. Below the critical pressure they rapidly but continuously change 
upon cooling. 

Computer simulations of the TIP5P and ST2 models show that response functions such 
as isobaric specific heat Cp and isothermal compressibility have maxima as functions of 
temperature if the system is cooled isobarically at P < P., Q. The loc, of these ma.- 
ima asymptotically approach one another [12| as the critical point is approached, because 
response functions become expressible in terms of the correlation length which diverges at 
the critical point. The locus of the correlation length maxima is called the Widom line. 

n 

Experimental studies of supercritical water [lq| indeed show that various response func- 
tions have sharp maxima in the region of the phase diagram above the liquid-vapor critical 
point Ci, but no direct experimental indication of a liquid-liquid critical point C2 had been 
available due to unavoidable crystallization of bulk water. It was found that water remains 
unfrozen in hydrophilic nanopores for T > 200 K |2^ . Moreover when cooled at constant 
pressure for P < Pc2 the dynamics changes from non-Arrhenius to Arrhenius at T = Tx (P). 
The line Tx (P) is located in the range of temperatures between 200 — 230 K and has a neg- 
ative slope in the P-T phase diagram. Computer simulations suggest that this line may be 
associated with Tvf(-P)) the Widom line, near which the local dynamic characteristics must 
rapidly change from those resembling the properties of HDL at high temperature to those 
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of LDL at low temperature [15|. 

Here we explore the hypothesis j21| that the observed glass transition in biomolecules 
is related to the liquid-liquid phase transition using molecular dynamics (MD) simulations. 
Specifically, we study the dynamic and thermodynamic behavior of lysozyme and DNA 



221 for fi) an 



24| 



in hydration TIP5P water, by means of the software package GROMACS 
orthorhombic form of hen egg-white lysozyme j^l and (ii) a Dickerson dodecamer DNA 
at constant pressure P = 1 atm, several constant temperatures T, and constant number of 
water molecules (NPT ensemble) in a simulation box with periodic boundary conditions. 
We first allow the system to equilibrate at constant temperature and pressure using the 
Berendsen method. This initial equilibration is followed by a long production run during 
which we calculate the dynamic and static properties. Equilibration times vary for different 
temperatures from a few ns for high temperatures to as much as 40 ns for low temperatures. 
The MD for DNA was performed using the Amber force field 25|. For lysozyme simulations, 
the system consists of a single protein in the native conformation solvated in = 1242 
TIP5P .olecu,. m Q. The. .n,.Mio„ condMon. _d to a .atio 

of water mass to protein mass of 1.56. The DNA system consists of a single DNA helix 
with 24 nucleotides solvated in = 1488 TIP5P water molecules, which corresponds to an 
experimental hydration level of 3.68. 

The simulation results for the mean square fiuctuations < > of protein are shown in 
Fig.^a). We calculate the mean square fiuctuations < > of protein from the equilibrated 
configurations, first for each atom over 1 ns, and then averaged over the total number of 
atoms in the protein. We find that < > [Fig.^a)] changes its functional form below Tj 



242 ± 10 K. Moreover, upon c oo. 
crossover from non-Arrhenius 



ing, the diffusivity of hydration water exhibits a dynamic 
to Arrhenius behavior at the crossover temperature Tx ~ 
245 ±10 K [Fig.^c)]. A similar temperature dependence of diffusivity of bulk TIP5P water 
was observed in Ref. The coincidence of Tx with Tp within the error bars indicates 

that the behavior of the protein is strongly coupled with the behavior of the surrounding 
solvent, in agreement with recent experiments [2l|. Note that Ty is much higher than the 
glass transition temperature estimated for TIP5P as Tg = 215K |lj|. Thus this crossover is 
not likely to be related to the glass transition in water. Here we will explore the possibility 
that instead it is related to a change in the properties of protein hydration water. 

We next calculate Cp by numerical differentiation of the total enthalpy of the system 
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(protein and water) by fitting the simulation data for enthalpy with a fifth order polyno- 
mial, and then taking the derivative with respect to T. Figure Efa) displays a maximum 
of Cp{T) at Tw ~ 250 ± 10 K for the case of lysozyme-water system j29|. The fact that 
Tp ^ Tx ~ Tw is evidence of the correlation between the changes in protein fluctuations 
[Fig. n(a)] and the hydration water thermodynamics [Fig. |21 (a)]. Thus our results are con- 
sistent with the possibility that the protein glass transition is related to the Widom line (and 
hence to the hypothesized liquid-liquid critical point). Crossing the Widom line corresponds 
to a continuous but rapid transition of the properties of water from those resembling the 
properties of a local HDL structure for T > Tw(-P) to those resembling the properties of a 
local LDL structure for T < Tw(-P) 20|. A consequence is the expectation that the fluc- 
tuations of the protein residues in predominantly LDL-like water (more ordered and more 
rigid) just below the Widom line should be smaller than the fluctuations in predominantly 
HDL-like water (less ordered and less rigid) just above the Widom line. 

To test this interpretation, we analyze the structure of hydration water on the two sides 
of the Widom line. Fig. Efa) shows the oxygen-oxygen radial distribution function g{r) on 
two sides of the Widom line for lysozyme hydration water. The first peak of g{r) on the low 
temperature (T=230 K) side is sharper and the first minimum is shallower compared to the 
g{r) on the high temperature (T=270 K, 300 K) side of the Widom line, suggesting that 
water is more structured on the low temperature side. Further, we calculate the structure 
factor S{q) of lysozyme hydration water [Fig. Etc)]. The first peak of the structure factor 
associated with the hydrogen bond is very sharp and pronounced, for T < Ty^{P), it is 
diminished and moves to larger wave vectors for T > Tw(-P), consistent with a LDL-like 
local structure for T < TyjiP) and a HDL-like local structure for T > T^^{P). 

Previous simulations Q| and experiments suggest a "glass-like" transition of DNA 
around temperature 230 K. Hence to test if the dynamic crossover depends on the solute, 
we performed a parallel study of the DNA Dickerson dodecamer [2^ . We find that fluctua- 
tions of the DNA molecule [Fig. ^b)] change their behavior approximately at the same 
temperature as lysozyme, with T ^ 247 ± 10 K. The dynamic crossover in the hydration wa- 
ter upon cooling from non-Arrhenius to Arrhenius behavior takes place at Tx ~ 250 ± 10 K 
[Fig. ^d)]. For DNA hydration water system Cp has a maximum at T ^ 250 ± 12 K, similar 



to the case of protein (see 
the DNA hydration water 



sys 

E 



ig-Hb)) |3l|. Fig.Efb) and Fig. Efd) show g{r) and S{q) for 
Further to describe the quantitative changes in structure 
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of hydration water, we calculate the local tetrahedral order parameter Q j33| for hydration 
water surrounding lysozyme and DNA. Figs. Efe) and El^f) show that the rate of increase of 
Q has a maximum at 245 ± 10 K for lysozyme and DNA hydration water respectively; the 
same temperatures where we find a crossover in the behavior of mean square fluctuations 
and a change in the behavior of the dynamics of hydration water. 

The quantitative agreement of the results for DNA and lysozyme suggests that it is 
indeed the changes in the properties of hydration water that are responsible for the changes 
in dynamics of the protein and DNA biomolecules. Our results are in qualitative agreement 
with recent experiments on hydrated protein and DNA j^| which found the crossover in 
side-chain fluctuations at Tp 225 K . 

Note Added in Proof: After this work was submitted, we learned of interesting parallel 
work on silicon, which also interprets structural change in g{r) and S{q) as crossing the 
Widom line 3^. 

We thank C. A. Angell, L. Cruz, P. G. Debenedetti, L. Liu, P. J. Rossky, and B. Urbane 
for helpful discussions of both experimental and computational work, and NSF for financial 
support. 
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FIG. 1: Mean square fluctuation of (a) lysozyme, and (b) DNA showing that there is a transition 
around Tp ~ 242 it 10 K for lysozyme and around Tp ~ 247 ± 10 K for DNA. For very low T one 
would expect a linear increase of < > with T, as a consequence of harmonic approximation for 
the motion of residues. At high T, the motion becomes non-harmonic and we fit the data by a 
polynomial. We determine the dynamic crossover temperature Tp from the crossing of the linear fit 
for low T and the polynomial fit for high T. We determine the error bars by changing the number 
of data points in the two fitting ranges. Diffusion constant of hydration water surrounding (c) 
lysozyme, and (d) DNA shows a dynamic transition from a power law behavior to an Arrhenius 
behavior at Tx ~ 245 it 10 K for lysozyme and Tx ~ 250 it 10 K for DNA, around the same 
temperatures where the behavior of < j;^ > has a crossover. 
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FIG. 2: The specific heat of the combined system (a) lysozyme and water, and (b) DNA and water, 
display maxima at 250 it 10 K and 250 it 12 K respectively, which are coincident within the error 
bars with the temperature Tp where the crossover in the behavior of < > is observed. 
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FIG. 3: (Color online) Oxygen-Oxygen pair correlation function g{r) for (a) lysozyme hydration 
water, and (b) DNA hydration water, on crossing the Widom line from the HDL side (T = 270 K, 
300 K) to the LDL side (T = 230 K). Structure factor of hydration water surrounding (c) lysozyme, 
and (d) DNA on two sides of the the Widom line. Upon crossing the Widom line, the local 
structure of water changes from more HDL-like to more LDL-like, reflected in the sharper and more 
prominent first peak. The first peak associated with the hydrogen bond distance also moves to small 
wave vectors, suggesting a change from the HDL to the LDL-like local structure of water at low 
temperatures. Derivative with respect to temperature of the local tetrahedral order parameter Q 
for (e) lysozyme and (f) DNA hydration water. A maximum in \dQ/dT\ at Widom line temperature 
suggests that the rate of change of local tetrahedrality of hydration water has a maximum at the 
Widom line. 



